M ensen zijn nieuwsgierig en wiiien graag ex¬ 
perimenteren. Soms voeren ze aiieriei experi¬ 
menten in huis uit met het gevaar van eiektrische 
schokken en chemische explosies. Tegenwoordig 
beschikken veien thuis over een computer waarmee ze, 
misschien zonder het nog te weten, een heei iaboratorium 
in huis hebben gehaald. Het is nu mogelijk vele ver¬ 
schijnselen van de natuur te begrijpen aan de hand van 


In dit boek beschrijft de auteur een groot aantal 
computerexperimenten waarmee men inzicht kan krijgen 
in zulke uiteenlopende zaken als het verloop van een 
epidemie, de groei van planten, de baan van een ruimte¬ 
vaartuig, milieuvervuiling en de werking van een 
pacemaker. Men kan zijn eigen planetarium maken, en 
men kan gaan begrijpen hoe bedrieglijk een Mercator- 
kaart kan zijn. Dat alles met eenvoudige formules, die 
niet meer dan een minimum aan wiskundige kennis 
veronderstellen. 

Alle experimenten zijn uitvoerbaar met eenvoudige 
computerprogramma’s, geschreven in de moderne en 
goed leesbare QBasic vorm. De stof is voor een deel 
ontleend aan de maandelijkse bijdragen van de auteur 
in het tijdschrift Natuur en Techniek. 
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1 Pythagoras en priemgetallen 


Over de historische figuur Pythagoras is maar weinig bekend. Het meeste 

een veelzijdig geleerde was, en dat hij zich in het bijzonder bezighield met 
de eigenschappen van de natuurlijke getallen, 1, 2, 3, 4, 5, ... De naar 
hem genoemde stelling over vierkanten beschreven op de zijden van een 
rechthoekige driehoek heeft hij kennelijk ontleend aan de wiskunde van 
de Mesopotamiërs. In bewaard gebleven kleitabletten in spijkerschrift 
vinden we daar bewijzen van. De stelling van Pythagoras is meetkundig 
geïllustreerd in figuur 1.1. 






























maar helaas werkt dit niet in sommige Basic-versies, omdat de mo 
opdracht slechts op gewone (korte) gehele getallen betrekking he 
Q(uick)Basic gaat het wel goed. In het programma PRIEM geb 
we de INT-opdracht: 


WHILE N\D = INT(N\D) ... 


Delen I 


ttijd 


Met de compilerversie Power Basic, de opvolger van Turbo Basic, is het 
zelfs mogelijk getallen bestaande uit achttien cijfers te factoriseren. Power 
Basic beschikt namelijk over de unieke mogelijkheid om intern met inte¬ 
gers van viervoudige lengte (quads) te werken. Het grootst toelaatbare 
gehele getal is dan 2*^-1, of in decimale notatie 

9 223 372 036 854 775 807 


DEFLNG J,F,T : DEFQUD N 

N een viervoudige lengte heeft. Bij Power Basic moet men we! bedacht 
zijn op een klein probleempje. Bij het inlezen van een geheel getal vanaf 
het toetsenbord zijn maar achttien decimale posities mogelijk. Het groot¬ 
ste getal dat in het hier geboden programma gefactoriseerd kan worden 
bestaat dus uit achttien negens, toch nog groot genoeg! 

Het is aardig eens een paar tijden te vergelijken. We nemen daartoe het 
eerder genoemde getal 2^'-l. Met de door de auteur gebruikte Compaq 
386/20 met coprocessor kost het in de moderne interpretatieve Quick- 

tien seconden voor nodig, terwijl zowel Turbo Basic als Power Basic er 

ziggehouden met de ontbinding van 100 895 598 169. Fermat vond de 
ontbinding 112303 * 898423. Met Power Basic gaat het nu in zes secon¬ 
den. Een recente computer van het i486/33 type presteert hetzelfde in 
slechts een enkele seconde, een resultaat dat de kracht van de moderne 
computer demonstreert maar dat ons ook met diep ontzag vervult voor 
het denk- en schrijfwerk van de twee wiskundigen, meer dan driehonderd 

lezers uit hun krachten te beproeven op de ontbinding van het getal 31 
059 639 803 977. Na 83 seconden meldt ons Power-Basic-programma 
het resultaat 


1597553* 19442009 




























2 Figuren van Lissajous 



















3 Wandelkaart en borduurpatroon 


Wie bij een wandeling in een berglandschap een goede wandelkaart ge¬ 
bruikt heeft veel plezier van de daarop aangegeven hoogtelijnen, niveau- 
lijnen of isohypsen. Langs een dergelijke lijn is gemeten ten opzichte van 
het zeeniveau de hoogte overal gelijk. De omgeving van een heuvel wordt 
bijvoorbeeld gemarkeerd door een aantal ellipsvormige niveaulijnen, 

1)! Heeft de top van de heuvel een vla4e vorm terwijl de hellin- 

liggen dicht tegen elkaar aan waar de helling steil is, en ze liggen ver uil 
elkaar wanneer we het vlakke topgedeelte naderen. Bevinden we ons op 
een bergpas, dan zouden we heel wat anders zien. Lopen we langs een 
pasweg en zijn we op het hoogste punt aangeland, dan zien we links en 

gaat. Op de kaart zien we dan een aantal hyperboolachtige bogen als geïl¬ 
lustreerd in figuur 3.1. 

Het is vrij eenvoudig hiervan met de computer een simulatie te maken, 
maar daarvoor is wel enige elementaire kennis van de analytische meet¬ 
kunde van de ruimte nodig. De hoofdzaak is dat we aan een functie/( xjj) 

helpt misschien wanneer we de variabelen x en y zien als beeldscherm- 
coördinaten van een Cartesisch coördinatenstelsel terwijl de a-as er lood¬ 
recht op staat, naar ons toegericht. In de analytische meetkunde van de 
ruimte beschikken we over een paar standaardoppervlakken als bol, ke¬ 
gel, paraboloïde, ellipsoïde en zadelvlak. Een bol met straal 1 is bepaald 







{xlay + {ylby + {z(cY^\ 
een (elliptische) paraboloïde door 
z = {xlay + (ylhy 

en een zadelvlak (hyperbolische paraboloïde) door 

z = (xlay-(y/by 

Voor de simulatie van een heuvel kunnen we een paraboloïde nemen, zeg 

is en dat de paraboloïde als een soort reflector van een fietslamp erboven 
uitsteekt. Om de bijbehorende hoogtekaart te maken gebruiken we het 
programma HOOGTE 1. 

In dat programma denken we het oppervlak uitgespreid over het vier¬ 
kant -I<x'<l,-l<3;<lals grondvlak. Dat vierkant verdelen we in een 
flink aantal punten, (2« + 1)^ voor w = 160 bijvoorbeeld, dus ongeveer 
100000. In elk van die punten bepalen we de hoogte z door z -y^ uit 
te rekenen. Het oppervlak verdelen we in horizontale plakjes, een aantal 















figuur 3.6 Een kruissteekjespatroon 


Men kan natuurlijk op goed geluk allerlei functies proberen, maar men 
kan ook de omgekeerde weg inslaan door aan te sturen op een bergland¬ 
schap dat een bepaalde meetkundige vorm dient te hebben. Men kan 
daartoe het landschap als het ware opbouwen uit bergen als elementaire 
bouwstenen. Een berg opgericht op het ^c^ij^-vlak met de top in de oor¬ 
sprong kan worden voorgesteld door de functie 

z = exp(-x*-y) 









^pC-ax^-bxy-cy^ 









4 Zwevingen en moiré 


in de frequenties heel we 


nig verschillen aanleiding geeft tot een periodiek aanzwellen en verzwak¬ 
ken van het geluid. Een pianostemmer maakt er handig gebruik van door 
bij het stemmen juist op dergelijke zwevingen te letten. Heeft de ene 
snaar bijvoorbeeld een frequentie van 500 trillingen per seconde en een 
naastliggende snaar er een van 502 Hertz dan constateert een goede luis¬ 
teraar, en een pianostemmer is zo iemand, tweemaal per seconde een ver¬ 
sterking en een verzwakking. Het is dus de kunst bij het stemmen zodanig 
te werk te gaan dat het zwevingseffect geheel verdwijnt. Het verschijnsel 
kan eenvoudig worden verklaard door twee hoogfrequente harmonische 
trillingen bepaald door de formules sin(ar) en sin(p) bij elkaar op te tel¬ 


len. De goniometrie leert ons dat 


sinCar) + sin(30 = 2 sinCCa + p)r/2)cosCCa - P)t/2) 


in hoofdzaak een harmonische trilling waarvan de frequentie het gemid¬ 
delde is van a en p, dat dus weinig van a en P verschilt. De amplitude 
2 cos((a-p)r/2) is daarentegen een langzaam veranderende, laagfrequen- 

mules graag bevestigd ziet door plaatjes op het computerscherm kan zich 
overtuigen met een computerexperiment. Daarbij kan men het program¬ 
ma ZWEVING gebruiken waarin de som van twee harmonische trillin¬ 
gen sin(27cs,r) en sin(27Ci20 grafisch wordt afgebeeld. Het resultaat ziet 





















































Een beetje theorie 

kennis die men in het vorige hoofdstuk nodig had om een hoogtekaart 
met niveaulijnen van een wiskundig oppervlak te kunnen begrijpen. Dat 


de hoogtelijnen eruit als concentrische ellipsen, en bij passen zijn het 
hyperbolische bogen. Is z de vergelijking van het oppervlak, dan 

zijn de niveaulijnen bepaald door = nh waarbij n de rij van de gehele 
getallen doorloopt en h meestal een klein getal is. Bij elke functie, en een 
keuze van /i, behoort een lijnenrooster. Het eenvoudige geval 

f(x^) = ax 

waarbij a een willekeurig getal is, geeft door gelijkstelling van ax met nh 
een stelsel equidistante evenwijdige lijnen bepaald door 







een stelsel equidistante concentrische cirkels. 

Het principe hoe roosterlijnen kunnen worden opgevat als niveaulijnen 
hebben we meetkundig geïllustreerd in figuur 4.6 voor een eendimensio¬ 
naal rooster waarbij de roosterafstand volgens een wortelformule ver- 







Het resultaat is dus 

Kxc!) = ihta^)^ 

hoe de parabool x als de doorsnede van een dal door equidistante 
horizontale niveaulijnen doorsneden wordt en hoe dat aanleiding geeft 


We beschouwen nog eens het eenvoudigste geval v; 
een frequentie a en een golflengte hia. We laten da 



l'(x)|. Het ligt voor de hand dat die zelfde uitspraak ook waar is wanneer 
/(x) en g{x) niet noodzakelijk lineair zijn. Voegen we de tweede variabele 
y er ook nog aan toe dan krijgen we de volgende fundamentele eigen- 


bepaald door |f(x,y) - gCxfy)|. 

We zijn nu in staat te begrijpen waarom we in het voorbeeld van MOI- 
REC2 als moiré-patroon een stelsel evenwijdige lijnen krijgen. De cirkels 




1 geeft 


Volgens de hoofdstelling is het patroon bepaald door het verschil 
f(.x+S,y)-f(.x,y)=25x+S‘^2Sx 

waarbij we een heel klein getal, verwaarloosd hebben. Gelijkstelling 
met nh geeft 
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cirkelL De programmals MOIRÉS 1 
RES4 mogen de lezer daarbij op weg 

























5 Potentialen en velden 





Een elektrische 









De tweedimensionale potentiaal in P(x,y) ten gevolge 



Ax^)=\og.{APIOF) 








f(x^) = -log(OP) = -logC;c^ +y)/2 

beeld die bij een enkele puntbron behoren. Het gebruikte programma is 
POTVELD1. 

met gelijke maar tegengestelde lading of intensiteit weergegeven, het be¬ 
kende patroon van twee elkaar loodrecht doorsnijdende cirkelbundels 
waarover we in de aanvang van dit hoofdstuk hebben gesproken. Het ge¬ 
bruikte programma is POTVELD2. 


wikkelder velden 





















die vraag als de vector van de partiële afgeleiden van ^ naar ;c en 3 ^. In 

met de plaatselijke daling in de ;i£:-richting en de overeenkomstige daling 
in dejv-richting 

dx =f(.x^) -f(x+k^), dj^ =fix^) ~f(x^+h) 

waarbij h een klein getalletje is. De grootheden dx en ds kunnen we, na 
het aanbrengen van een evenredigheidsconstante, interpreteren als de 
verplaatsing van het pimt dat zich langs een stroomlijn over het beeld¬ 
scherm beweegt. In figuur 5.7 laten we zien hoe het in figuur 5.6 afge- 


6 Zelfgelijkvormige plantengroei 


Biologen hebben al vaak opgemerkt dat de bouw van levende wezens niet 
alleen veel regelmatigheid vertoont, maar dat het soms is alsof een zelfde 
groeiproces zich op verschillende plaatsen en tijden telkens herhaalt. We 
hoeven maar naar een boom, een vogelveer of een varenblad te kijken om 
daarvan overtuigd te raken. Een van de grondgedachten van de natuur¬ 
wetenschap is om te trachten de veelsoortige verschijnselen op 20 een¬ 
voudig mogelijke wijze te kunnen begrijpen. We plaatsen ons als het ware 
op de stoel van de Schepper en we gaan uit van een paar simpele natuur¬ 
wetten van meetkundige aard. We zijn daarbij in het goede gezelschap 
van de Griekse filosoof Plato, die reeds verklaarde dat God bij de schep¬ 
ping van de wereld op meetkundige wijze te werk was gegaan. Om een 
idee van de mogelijkheden te krijgen en in een streven om alles zo een¬ 
voudig mogelijk te houden beperken we ons in dit hoofdstuk tot een we¬ 


reld va: 


it begrip van de gelijkvormig- 

itsing en een inkrimping. Een 
;en gevolge van een dergelijke 
is gedraaid en vervormd, te- 


te beginnen met een vel papier waarop bijvoorbeeld een vierkant is gete¬ 
kend. Daarop leggen we een transparant vel, en we kopiëren daarop het 
vierkant. Het doorzichtige vel verschuiven we en tegelijk draaien we het 
wat. We kunnen nu goed het verplaatste vierkant vergelijken met het ori¬ 
gineel in het onderste vel. Het bovenste vel kunnen we ook in zijn geheel 
omdraaien en weer in een andere positie op het onderste papier leggen. 
We hebben dan te maken met een spiegeling. Voor een vierkant maakt 
dat niets uit, omdat een vierkant zelf spiegelsymmetrisch is, maar zouden 
we als basisfiguur bijvoorbeeld de letter F hebben genomen, dan is er na¬ 
tuurlijk wel verschil. Om het effect van een inkrimping te zien kunnen we 
het beste beginnen met het aanbrengen v 

speciaal, affiene transformaties, 
altijd recht blijven, en dat een w 


C6.1) 


cx^dy+f 




voorbeeld een wülekeurige rechthoek. Het oppervlak van de rechthoek 

langrijk getal. Is dat getal negatief dan heeft de transformatie het karakter 
van een spiegeling, is dat getal positief dan is er een gewone vervorming. 
Is ten slotte ad-bc = Q dan is de transformatie een projectie waarbij alles 
tot een enkele rechte lijn wordt samengedrukt. 

Bij een aantal concrete voorbeelden heeft men te maken met een paar 
speciale gevallen van (6.1) die elders uitvoerig besproken zijn, maar waar¬ 
over we hier toch iets moeten zeggen. Wanneer bij een gelijkvormig- 


ot de oorsprong verkort 


it. Zo’n pi 


De oorsprong (0,0) blijft op zijn plaats en is dus een dekpunt. Is in het 
bijzonder a = 0 dan reduceert (6.2) zich tot 


Ie verkortingsfactor van lijnstukken, wel te onderscheiden van de tweedi¬ 
mensionale verkortingsfactor van rechthoeken en dergelijke die gelijk is 


De boom van Pythagoras 

lijkvormige plantengroei is de beroemde boom van Pythagoras, ontstaan 
in de Tweede Wereldoorlog op de tekentafel van ir. A. Bosman en afge- 





























i de Pythagoras-boom, het resultaat van twee contracties L en R, gaat 
hele figuur op gelijkvormige wijze over in een deel van de figuur. 



In figuur 6.8 zien we hoe het varenblad als een spel van chaos en orde 
kan worden samengesteld uit vier gelijkvormigheidstransformaties. We 
gaan daarbij uit van een smalle rechthoek geplaatst aan het begin van het 












7 Wiskunde en wiskunst 
















Kunst op de wijze van Mondriaan 

We beginnen met het programma MONDRIAN waarin nog geen kleu- 

voorgeschreven aantal horizontale of verticale lijnstukjes getekend waar¬ 
bij plaats en grootte volgens toeval bepaald zijn. De opdracht RANDO- 
MIZE #, waarbij voor # een willekeurig geheel getal gekozen kan worden. 























de waarde te geven. De grootte van het gekleurde vierkantje kan met de 
gekleurd blokje in het beschikbare veld. Hel aantal kleuren hangt af van 
0.7 betekent dat ongeveer zeventig procent van de velden gekleurd 

Wanneer op het scherm het patroon verschijnt kan met een willekeuri¬ 
ge toetsaanslag een andere kleuring worden verkregen. Dat kan eindeloos 
worden herhaald, en slechts met de ‘X’-toets kan het programma worden 

‘othercolours’, in principe het gebruik van de palette-opdracht. Men zij 
er evenwel op bedacht dat Turbo Basic en Q(uick)Basic zich ten aanzien 
van de palette-opdracht in de VGA-modus verschillend gedragen. Zie 
figuur 7.2, bladzijde 83 in het kleurkatern. 


Kunst op de manier van Klee 

stellingen gemaakt waarbij bepaalde elementaire figuren op schijnbaar 

gaafde kinderen zijn geweest. In het programma KLEE betuigen wij 

84). ^ ^ 

Het programma is verwant aan MONDRIAN, maar in plaats van ho¬ 
rizontale en verticale lijnstukjes worden nu op willekeurige plaatsen ge¬ 
kleurde lijnstukjes of driehoekjes in een willekeurige oriëntatie en van een 
door toeval bepaalde grootte geplaatst. Het programma is zodanig ge- 

Men zou bijvoorbeeld de subroutine ‘lijn’ kunnen vervangen door een 
subroutine ‘cirkel’ met behulp waarvan gekleurde cirkelschijfjes kunnen 
worden gemaakt, bijvoorbeeld 

C=1+INTC4*RND) 

CIRCLE (X,Y),V,C : PAINT (X,Y),C 
RETURN 


Oude volksgebruiken dreigen te verdwijnen omdat onze moderne maat¬ 
schappij zo sterk op consumptie van in overvloed aanwezige massawaar is 
gericht. Toegegeven moet worden dat maar weinigen over artistieke ver- 

paaseieren te beschilderen. In elk geval zijn er vele mogelijkheden om op 
het beeldscherm patroontjes te ontwerpen die kunnen dienen als voor¬ 
beeld van beschildering van een paasei. Het eerste probleem dat zich aan- 
60 dient is de wiskundige simulatie van een ei. 


Een ei is wiskundig gesproken een omwentelingslichaam. Doorsneden 
loodrecht op de omwentelingsas zijn altijd cirkels. Een doorsnede die de 

tijd eender uit. Het is voldoende wanneer we zo’n doorsnede bekijken. 
We gebruiken een x,3;-coördinatenstelseI waarbij we de x-as als omwente¬ 
lingsas kiezen. Die as is dus een symmetrieas van de figuur. Het eenvou¬ 
digste ei is een zuivere bol, en ook om biologische redenen is dat de beste 
vorm, omdat een bol bij een gegeven inhoud het kleinste oppervlak heeft. 
Maar omdat een ei in de natuur door een cilindrische buis moet worden 
getransporteerd zijn eieren van nature enigszins cilindrisch samenge- 

De cirkelvorm is wiskundig bepaald door 

y = of y=±^(l-x^) 

dig te verkleinen, y-^ay waarbij a < 1. We krijgen dan de functionele be- 



y=±a'l(l-x^)ll + bx) 

een bevredigend resultaat geeft. In figuur 7.4 zien we dat de eivorm heel 
goed benaderd wordt. De simulaties kunnen met het programma EI wor¬ 
den uitgevoerd door in het programma de waarden van a en van b telkens 









8 Het pad van de schildpad 



het nooit lukken. 

Leonhard Euler (1707-1783) was in zijn tijd een der grootste en meest 
produktieve wiskundigen. Zowel zuivere als toegepaste wiskunde hadden 
zijn aandacht. In zijn Lettres d une princesse d’Allemand toonde hij ook de 
wetenschap op vaardige wijze te kunnen populariseren. In zijn tijd 
stroomde in het toenmalige Koningsbergen de rivier de Pregel waarin 
zich tussen beide oevers een eiland bevond. De stadsdelen en het eiland 
waren onderling verbonden door zeven bruggen (figuur 8.2). Men vroeg 

dat elke brug slechts eenmaal gepasseerd werd en dat men weer tot het 
punt van uitgang zou terugkeren. Euler liet zien dat dat niet mogelijk was. 
Hij simuleerde de wandeling door een graaf met knopen en trok uit de 

werd hij zo de schepper van de grafentheorie. 


De schrijvende schildpad 

Ouder geworden beschikken we over modem tekengereedschap en dat is 
tegenwoordig de computer, een machtig hulpmiddel dat veel mogelijk 

wennen, en een van de betere methoden daarbij is het gebruiken van wat 
in de computerwetenschap turtlegraphics heet. 

De kinderen zien dan een soort schrijvende schildpad zich over het 
beeldscherm bewegen die met eenvoudige toetsinstructies te besturen is: 












len als F, + en De schildpad kenmerken we door een drietal getallen 
X, y en <p. De getallen ;c en 3^ geven de positie van de schildpad in recht¬ 
hoekige coördinaten waarvan we de oorsprong vaak in het centrum van 

opzichte van een verticale of horizontale lijn door het centrum. Het derde 
getal geeft de richting waarin de schildpad zich beweegt. Werken we in 
graden dan correspondeert (p = 0 met horizontaal naar rechts, (p = 90 met 
verticaal omhoog enzovoort. De opdrachten + en - betekenen dan dat (p 
telkens met een bedrag van 90“ vermeerderd of verminderd wordt. 
Eventuele grote waarden van de hoek kunnen we natuurlijk altijd redu¬ 
ceren tot een waarde tussen zeg 0 en 360'. Een hoek van bijvoorbeeld (p = 

We spreken af dat we aan het begin de toestand (jc^yatp) van de schild¬ 
pad aangeven als (0,0,0). De opdrachten + en - kunnen we wiskundig be¬ 
schrijven als 



lengte in de richting die de schildpad dan heeft. Die lengte moet van ge¬ 
in het programma TURTLE is het een en ander nader uitgewerkt. 
Het programma bestaat eigenlijk uit een telkens te herhalen lus die re¬ 
ageert op een toetsaanslag. Met de toets F schrijft de schildpad een recht 
lijntje. Met de plus-toets wordt de richting veranderd in een kwartslag 
naar links (in positieve zin volgens de wiskundige conventie) en met de 
min-toets draait de schildpad een kwartslag naar rechts. Aldus kunnen 
we met een aantal toetsaanslagen de schildpad een uit horizontale en ver¬ 
ticale lijntjes opgebouwde weg laten afleggen. Het programma laat zien 
dat met een X- of x-aanslag het einde is bereikt. Heel vaak worden in een 
programma de toetsen X, x, Q, q gebruikt om te eindigen (X = eXit, Q = 
Quit). De staplengte is aan het begin als 10 vastgesteld, maar dat getal 
kan natuurlijk veranderd worden. Gebruikers van een oudere Basic-ver- 

programmeertaal, een betere lusstructuur, een handige meerkeuzetoets 


De taal van de schildpad 

We hebben waarschijnlijk al gemerkt dat de schildpad best tot een reeds 
afgelegd stukje weg kan terugkeren. Bezwaarlijk is dat niet, maar mis- 




De richting van de schildpad wordt dan nog steeds door de hoek (p ge¬ 
markeerd, maar in plaats van de hoek cp gebruiken we het tweetal (p,g), de 











opgave voor de lezer om uit te rekenen hoe groot de volgende woorden 
worden. Maar het is inmiddels wel duidelijk dat we de computer hard no- 



graphics zien. De techniek van turtlegraphics is in principe heel een- 

een lang woord, een string, gevormd dat alle instructies bevat die nodig 
zijn om de schildpad de verlangde figuur te laten vormen. In de praktijk 

































































F-F (ofF + F) 


geeft precies hetzelfde resultaat, zij het dat de lijn van Lévy in zijn geheel 
om het beginpunt geroteerd is. Is p de gebruikte orde, dan is de totale 
draaiboek gelijk aan p * 45°. Het grote voordeel is dat de string veel min¬ 
der snel aangroeit zodat p= 12 normaal gesproken gemakkelijk haalbaar 

De tweede vorm van de winkelhaak, geïllustreerd in figuur 8.11b, de 
vorm die ten grondslag ligt aan de drakenlijn, geeft aanleiding tot twee 
produktieregels 


G->-F + + G- 


ZX pzi 
I7Z 



.F-G 
► F + G 


waarbij voor + en - ee 
Het met deze rege! 


gramma TURTLED kan dan snel 
ï groeifase, op het scherm toveren. 












9 Lijnen die een vlak vullen 


gemakkelijke stoel, naar televisiebeelden kijken beseffen we niet of nau¬ 
welijks wat er allemaal komt kijken voor de overdracht van een enkel 
beeld. Een beeld is begonnen als een tweedimensionaal plaatje, wordt ge¬ 
codeerd als een eendimensionaal signaal dat via antenne of kabel in onze 
apparatuur wordt opgevangen, en daarin wordt op een omgekeerde wijze 
het oorspronkelijke beeld weer samengesteld. ïn abstracte termen gaat 

we de punten van een rechthoek of een vierkant afbeelden op de punten 





drakenlijn. De door Peano gevonden lijn is natuurlijk een wiskundige 
idealisatie, maar een benadering van een dergelijke lijn is in de praktijk 
heel goed bruikbaar. Na Peano hebben andere wiskundigen soortgelijke 

wiskundigen bij de aanvang van de twintigste eeuw. Gecombineerd met 
de moderne computertechniek kan het gebruik van een Peano-lijn leiden 
tot een doeltreflfende methode om op systematische wijze een fysisch op- 

Afgezien van de theoretische en technische merites van Peano-achtige 
lijnen zijn ze een dankbaar object voor de computer. We zullen dan ook 
laten zien hoe de in het vorige hoofdstuk geïntroduceerde schildpad op 
grond van een paar eenvoudige regels een lijn van Peano kan beschrijven. 
Als voorbeeld nemen we de versie van David Hilbert (figuur 9.1). 

De schrijvende schildpad kan met behulp van het axioma X en de twee 
produktieregels 

X^-YF + XPX + FY- 
Y-^ + XF-YFY-FX + 

bool F is het tekenen van een lijnstuk, + betekent een kwartdraai naar 
links (of rechts) en - een kwartdraai naar rechts (-of links). De symbolen X 

de schildpad bij het tekenen niet gebruikt. De lezer doet er verstandig aan 
om zelf, gewapend met potlood en ruitjespapier, een klein stukje mee te 
tekenen. Men begint met het potlood ergens in een roosterpunt te plaat- 

ging. De eerste fase levert volgens de eerste produktieregel de string 
-YF+XFX+FY-. Staat de schildpad bij de startpositie naar rechts gericht, 

weer omhoog. Ten slotte heeft de schildpad weer de oorspronkelijke 
oriëntatie. In de tweede fase worden de symbolen X en Y volgens de 
produktieregels elk vervangen door een deelstring van elf symbolen. De 
totale string bestaat nu uit 51 symbolen. Daarin zitten vijftien symbolen 
F, zodat het meetkundige equivalent een gebroken lijn is van vijftien 
stukjes. De symbolen X en Y hebben geen directe meetkundige beteke- 
































































Regdl X -» XFYFX+F+YFXFY-F-XFYFX 
Regel2 Y -» YFXFY-F-XFYFX+F+YFXFY 








Figuur 7.2 Een toevallig geblokt dessin 













Figuur 15.9 c,a-Diagram van fasekoppeling 





F[-F][+F[+F]F]F[-F][+F] 



een hoek (p, die de looprichting aangeeft. Zodra de schildpad een vertak¬ 
king, een knoop van waaruit verschillende voortzettingen mogelijk zijn, 
tegenkomt dient hij dus die drie getallen te noteren. Hij doet dat op een 
soort kaartje - denk maar aan een speelkaart. Telkens wanneer de schild¬ 
en legt hij het nieuwe kaartje bovenop het oude. Op die wijze kan een sta¬ 
peltje van kaarten ontstaan, iets wat in computerjargon als een stack be¬ 
schreven wordt en wat wij gewoon een stapel noemen. Zodra de schildpad 
met het symbool ] het einde van een tak bereikt behoeft hij slechts naar de 
bovenste kaart van de stapel te kijken om te weten waar hij moet terugke¬ 
ren en welke richting hij dan moet inslaan. De kaart is dan niet meer no¬ 
dig en wordt verwijderd. We zien aldus in onze verbeelding een stapel 
kaartjes van wisselende samenstelling ontstaan. Het enige wat er gebeurt 
is een kaart op de stapel leggen of een kaart van de stapel weghalen. De 
stapeltechniek wordt vooral gebruikt bij de methode van backtracking, 
die in vorige hoofdstukken een paar maal ter sprake is gekomen. In de 
techniek van het programmeren wordt een stapel veelvuldig gebruikt. 
Het hart van menige computer, de centrale microprocessor zoals die van 


voeren. Moeilijk is dat nietj en het programma LIND1 laat zien hoe dat in 
zijn werk gaat. Eerst zetten we de verschillende symbolen van de schild- 
padtaal nog even op een rijtje. Met de positie van de schildpad bedoelen 
we daarbij de twee plaatscoördinaten x^y&n de bewegingsrichting (p. 












besproken model eraan toegevoegd. 


Axioma: X (figuur 10.4) 

Regel 1: X F [+X] F [-X] +X 
Regel2rF^FF 


Rege 

Rege 


a:Y (figuur 10.5) 

: Y ->YFX[+Y][-Y] 
:X-^X[-FFF][+FFF]FX 


Axioma: X (figuur 10.6) 

Regel 1: X ^ F- [[X] +X] +F [+FX] -X 
Regel2: F -> FF 


Zoals we in de vorige hoofdstukken hebben opgemerkt krijgen we proble- 
laten groeien. De string waaruit de opbouw van het groeisel bij de grafi- 

schikbare stringgeheugen. Maar ook nu weer kan het veel handiger door 
gebruik te maken van de techniek van backtracking. We krijgen in het 
resulterende programma nu twee soorten stapels. De eerste stapel dient 
alleen voor de opbouw van de string in de hoogste groeifase, de tweede 
stapel heeft alleen betrekking op de grafische bewerking om de boekhou¬ 
ding van de vertakkingen te verzorgen. 

Het programma LIND2 is aldus een versie van het programma 
LIND1 waarbij nu backtracking is toegepast. Het heeft in de in het boek 

andering van axioma en regels en met aanpassing van de beeldbegrenzing 


iurlO.3 Een plant met de regels 
FF+l+F-F-Fl-ï-F+F+F] 















11 Kartografie, het bedrog 
van landkaarten 

















waarbij a met de meridiaanlengte en P met de breedte correspondeert. 
Voor P = 0 krijgen we punten op de evenaar, voor P = 90° de noordpool 

Allereerst stellen we ons de opgave op het beeldscherm een afbeelding 
van de bol met een aantal breedtecirkels en meridiaancirkels te verkrij¬ 
gen. De eenvoudigste methode is die waarbij punten (x^,z) van de bol 
loodrecht op het beeldscherm geprojecteerd worden. Daarbij voeren we 















lipsvormige stippellijn. Aan het begin van het programma hebben we de 
hellingshoek (p van de rotatieas de willekeurige waarde van 25° gegeven, 
maar het is leuk hetzelfde programma te herhalen met een andere waarde 
van (p. Overigens wordt in het programma met een lusstru 
tal breedtecirkels getekend. Op overeenkomstige wijze w( 
meridiaancirkels gevormd. Het resultaat is weergegeven in 


fis 



Cilinderproje 










i I i' 


; = logtan(7ï/4 + p/2) 


Mercator zou deze formule niet hebben kunnen gebruiken, omdat de lo- 



conform, maar niet equivalent. Om een indruk van de grootte van be¬ 
paalde gebieden te krijgen is de Mercator-kaart nagenoeg onbruikbaar; 
een gebied als Groenland lijkt bedrieglijk groot als ware het een continent 










is. Het land denken we gelegen tussen 50° en 70°, maar dat betekent niet 
dat zuid- en noordgrens ook de bijbehorende breedtecirkels zijn. De wer¬ 
kelijke grenzen zijn bogen van grote cirkels en die worden bij de kaartpro- 


12 Tussen de sterren 



(12.1) 


m.d^r/dï^' =-Gm^mJir^ - 




waarin de gebruikte symbolen de volgende betekenis hebben: r, en zijn 

Voor de tweede ster geldt een overeenkomstige vergelijking, maar 
nieuwe informatie geeft dat niet omdat het zwaartepunt van de twee 
puntsterren steeds met de oorsprong samenvalt, in formulevorm 



(12.4) 


r^=aJ(l+ecos(p) 










Voor (p = 0 is de voerstraal zo klein mogelijk. De planeet staat dan zo 
dicht mogelijk bij de zon en de afstand bedraagt a/(l+e). Voor (p = 7t heeft 
de planeerde grootste afstand tot de zon en die is dus al{\-e). De som van 
beide grootheden, namelijk is de grootste diameter van de pla- 

netenbaan. De tweede vergelijking 

rMcp/dt = constant 

geeft de mogelijkheid om aan de baan een tijdschaal te verbinden. In het 
programma PLANEET geven we een aantal planetenbanen die op deze 
twee vergelijkingen gebaseerd zijn. Er worden daarin een achttal plane¬ 
tenbanen berekend en afgebeeld. De excentriciteit neemt telkens een 
beetje af. We beginnen met een sterk ellipsvormige baan en eindigen met 
een baan die bijna een cirkel is (figuur 12.2). 


meer verre planeten te verkennen. Een onderneming waaraan behalve 
technisch vernuft ook heel nauwkeurig rekenwerk ten grondslag dient te 

kracht van de aarde te laten overwinnen, en er is maar weinig energie over 
om onderweg wat bij te sturen. We lezen dan dat de baan zo gekozen is 
dat de satelliet na het verlaten van de aarde door een planeet op dusdani¬ 
ge wijze wordt aangetrokken dat die er als het ware wat extra energie aan 
toevoegt, zoals een kogelslingeraar met een extra draai van zijn lichaam 
een kogel ver weg weet te slingeren. 



















13 Maak uw eigen planetarium 






s(<p) cosca). 


sin(<p) cos(a) 


De constante c is de hoeksnelheid waarvan de waarde bepaald is door de 

waarnemen is het oneindig verre punt van het lijnstuk diens richting 
dus, waarbij P de aarde en P, de planeet is. Het komt er dus op neer op de 
een of andere manier die richting als punt op het beeldscherm weer te ge¬ 
ven. In het programma PLANEET2 is het een en ander gerealiseerd. De 
banen van aarde en Jupiter zijn zichtbaar gemaakt als cirkels in het 
vlak waarbij nu de inclinatie van de baan van Jupiter is verwaarloosd. De 
vector PP, wordt berekend en vervangen door een eenheidsvector, omdat 

de bovenkant van het beeldscherm de projectie op het x,2-vlak als punt 
afgebeeld. Het is dus alsof we in het vlak van de aardbaan naar die vector 
kijken. Door de beweging van aarde en Jupiter beweegt die projectie heen 
en weer. Zou de baan van Jupiter met die van de aarde samenvallen, dan 
zou die heen en weer gaande beweging volgens een rechte lijn gebeuren, 
niet zo interessant. Maar juist omdat de baan van Jupiter een kleine incli- 






natie heeft, beweegt het beeldpunt zich ook een beetje naar boven en naar 
hemelbeschouwer waarneemt bij het observeren van de planeten. Overi- 

het enige doel het verkrijgen van een wat beter plaatje is. De totale tijds¬ 
duur is heel willekeurig op vijftig radialen, dus zowat acht jaar, gesteld 
(zie figuur 13.2). 






14 Relaxatietrillingen 


Trillingsverschijnselen in de natuurwetenschappen kunnen op verschil¬ 
lende wijzen gemodelleerd worden. Het eenvoudigste model is de harmo- 

schreven. We kunnen u opvatten als de uitwijking van een deeltje uit de 

beeld. De constante ct is de frequentiej waarbij 2ït/o de trillingstijd is, de 
tijd die nodig is om precies een enkele trillingsgolf uit te voeren. Het aan¬ 
tal trillingen per seconde is dus a/27ü. 



(14.2) 


nentieel af volgens de factor expC-|ir/2). 



Beginnen we met een gegeven startpositie u=u^tnv= v^, in fysische ter¬ 
men een gegeven uitwijking aan het begin en een gegeven beginsnelheid, 
dan kunnen we op de gebruikelijke wijze met de opvolgers een baan 
vormen. De index n heeft nu de betekenis van het tijdstip nh gerekend 
vanaf het begin. In de bi jbehorende grafische voorstelling hebben we ken- 



























In het programma POL2 worden positie u en snelheid v tegen elkaar uit¬ 
gezet. Wanneer = 0, dan is de vergelijking van Van der Pol identiek met 
die van de ongedempte harmonische slinger. In plaats van (14.8) krijgen 









15 Het hart als elektrisch systeem 










































m(l^)op 




in (x,0)va 


sn (0^).Hetz 


t ruimtelijk inzicht heeft kan het volgende gedachtenexpe¬ 
riment uitvoeren. Hij neemt als hoofdvierkant een vierkant vel papier en 

ker. Boven- en onderrand van het vierkant worden daarbij aan elkaar 
vastgeplakt. De koker wordt nu in horizontale zin rondgebogen - het pa¬ 
pier rekt of krimpt daarbij wat mee - en de twee cirkelvormige uiteinden 


ject dat doet denken aan een werpring of een binnenband. Amerikanen 
die meer vertrouwd zijn met produkten als hamburgers en doughnuts 
worden altijd getroffen door de gelijkenis van het meetkundige lichaam 
met hun populaire lekkernij. Overigens schijnt ook hier de ‘donut’ inge- 
burgerd te raken. We zouden dan aan figuur 15.5 nog een grafische voor¬ 
stelling kunnen toevoegen waarbij de rechte lijn zich als een kronkellijn 
over het torusoppervlak beweegt. 

Uiteraard is tekenen op een plat vierkant heel wat eenvoudiger dan te¬ 


kenen op een torus, gevolgd door een projectie op een plat vlak. Maar 
toch kan dat torustekenen, denkbeeldig of desnoods schetsmatig, ons 
helpen inzicht te krijgen in fasemodellen. 


Torustekenen 

Het programma TORUS stelt de lezer in staat om op het beeldscherm 
de schijnbare omtrek van een torusoppervlak te verkrijgen en daarover 
een willekeurige ‘rechte’ lijn te laten kronkelen. Van de wiskundige ach¬ 
tergrond vermelden we aüeen dat het torusoppervlak in ruimtelijke coör- 





















i waarbij 


De coördinaten van een punt op de torus zijn gegeven door t en u 
t/(2K) en u/(2k) de betekenis van fase hebben en daarmee overee 
men met de horizontale coördinaat x en de verticale coördinaat van het 
vlakke periodieke vierkant. De rechte lijn van figuur 15.5 zou dan op de 
torus door u = 0.75t bepaald kunnen zijn. 

Het gaat eigenlijk nog steeds om een grafische illustratie van het feno¬ 
meen van Wenckebach, de fasekoppeling van 4 : 3, waarbij de horizonta¬ 
le coördinaat x van het periodieke vierkant en t van het torusoppervlak 
viermaal rondgaat terwijl de verticale coördinaat y of u drie keer rond¬ 
loopt. In plaats van 4 : 3 kunnen we elke andere verhouding nemen^ bij¬ 
voorbeeld p : q waarbij p on q gehele getallen zijn en piq een onvereen- 
voudigbare breuk is. Op het torusoppervlak komt dan een gesloten kron¬ 
kellijn die horizontaal ^maal en verticaal g-maal rondloopt. Wanneer de 

aan chaos oproept. In figuur 15.6 is om illustratieve redenen gekozen 
voor een fasekoppeling van 1:11. 


In 1914 was een jong fysioloog aan de McGill University School of 
Medicine van Montreal bezig met experimenten waarbij hij onderzocht 
hoe een han reageert op kleine elektrische prikkels. Hij had een apparaat 
gebouwd waarmee hij kleine stroomstootjes kon opwekken met een 
nauwkeung regelbare sterkte en frequentie. Experimenten met proefdie- 
















is die van de fase aan het begin en aan het eind van de cyclus. Het model 
is een fraai voorbeeld van een fase-fase-afbeelding die zich goed leent 
voor afbeelding op een torusoppervlak of, eenvoudiger, op een periodiek 


Simulaties 


Met het pacemaker-model kunnen we verschillende computerexperi- 
menten uitvoeren. Daartoe is een aantal programma’s beschikbaar met 
de gemeenschappelijke titel CIRKELH. 

In het eerste programma, CIRKELH 1, wordt het gedrag van het 

Eenvoudig is te zien dat er in sommige gevallen een aantrekkend dekpunt 
of een periodieke aantrekkende cyclus bestaat, fasekoppeling dus. In vele 
gevallen is er chaos, maar dat is in dit model niet goed waar te nemen. 

In het tweede programma, CIRKELH2, wordt de fase telkens tot 
(0,1) gereduceerd en wordt in het iteratieproces alleen het hoofdvierkant 
weergegeven. Gevallen van fasekoppeling en van chaos zijn hierin duide- 

In het programma CIRKELH4 kiezen we wederom a vast (bijvoor- 



hier merken we slechts op dat het patroon ten opzichte van c = 0.5 sym¬ 
metrisch is. Het diagram behorende bij de waarde a = 1 is afgebeeld in 

Het programma CIRKELH5 geeft een fraai diagram van fasekoppe¬ 
ling waarin gebruik van kleuren wel noodzakelijk is. De rechthoek be¬ 
paald door de grenzen 0<c<len0<a<27i correspondeert op het 
beeldscherm met een rechthoek van pixels. Voor elke pixel, equivalent 
aan een punt (c,a), kan een periode 9, de noemer van een eventuele fase¬ 
koppeling p : q, worden berekend. Vaak is q een niet te groot getal, maar 
wanneer dat getal hoog is en er dan een mogelijkheid van chaos bestaat, 
wordt q = 0 genomen. Volgens de kleurencode l+(? modi 4) wordt hier¬ 
uit een kleur van de pixel afgeleid. Wanneer het programma klaar is, na 
vele uren wellicht, is een fraai diagram verkregen dat is afgebeeld in 
figuur 15.9 (zie bladzijde 87 in het kleurkatern). 































16 Orde en chaos 


















De wiskundige beschrijving is 


het hoofdvierkant met de hoekpunten C0,0),ClsO),Cl,l),CO,l) over in het 
parallellogram met de hoekpunten C0,0),(l,l)s(2,3),Cl,2). Dat parallel¬ 
logram heeft dezelfde oppervlakte als het oorspronkelijke vierkant. Maar 








eigenlijk moeten we de transformatie van Amold opvatten als een trans¬ 
formatie op een torusoppervlak, zoals we die in het vorige hoofdstuk 
hebben besproken. Het hoofdvierkant staat dan model voor het gehele 
torusoppervlak. De mod 1 toevoeging betekent dat alles wat door de 
transformatie buiten het hoofdvierkant terechtkomt daarin weer wordt 
teruggebracht. We werken in feite op periodiek ruitjespapier als bespro¬ 
ken in het vorige hoofdstuk. 

Arnold laat zien hoe - zoals in figuur 16.2 - de kop van een kat, opge¬ 
bouwd uit een paar eenvoudige meetkundige lijntjes, vervormd en ver¬ 
knipt wordt. Oorspronkelijk bevindt de kattekop zich in het hoofdvier¬ 
kant. Door de transformatie gaat het hoofdvierkant over in een parallello¬ 
gram en komen delen van de kat in andere vierkanten terecht. Die delen 
moeten naar het hoofdvierkant worden teruggebracht. Dat betekent 
meetkundig dat het parallellogram in vier stukken wordt .geknipt en dat 
de vier delen als een legpuzzel op de plaats van het hoofdvierkant aan el¬ 
kaar worden gevoegd. Uiteraard is daarbij de kattekop in vreemde stuk¬ 
ken uiteengevallen. 

Wanneer de kat van Arnold op iteratieve wijze enige malen behandeld 

zijn geraakt. Dat is ook de bedoeling van Amold: laten zien hoe bij een 
heel eenvoudig model reeds chaos en volledige menging kan ontstaan. 
Wanneer de transformatie van Amold maar vaak genoeg wordt herhaald 
bevat elk stukje van het hoofdvierkant, hoe klein ook, wel onderdeeltjes 
van de versnipperde kat. Dat zelfde geldt natuurlijk ook bij de interpreta¬ 
tie van de kat op het torusoppervlak. 

Zowel bij de bakkerstransformatie als bij de ‘cat map’ is het dynamisch 
systeem chaotisch in de volgende betekenis. Twee banen die tegelijk be- 




lit elkaar en gaan steeds meer C' 
lar te maken willen hebben, e< 
uitdmkkingswijze die desgewenst vervangen mag worden door de const 
tering dat de onderlinge correlatie verdwijnt. 

Arnold vermeldt ook een eenvoudig model waarbij wel sprake is v; 
een goede menging, maar niet van echte chaos. Het betreft een torusi 




;n liggen op de kronkellijn 

























































Het overeenkomstige programma LYAPUN2, dat maar heel weinig 
van het vorige programma verschilt, geeft de grafiek van de Lyapunov-ge- 
tallen. Een kleine bijzonderheid is dat de waarde nul bij de gemiddelde 

nov-getal van min oneindig. In de grafiek zien we hooguit een plaatselijke 

Wanneer het systeem ordelijk is kan de gemiddelde groei soms op een¬ 
voudige wijze vooraf bepaald worden. Voor 1 < a < 3 is er een aantrek¬ 
kend dekpunt waar bijna alle banen naar toegaan. Dat dekpunt is 1-1/a 
en de groeifactor voor dat dekpunt is volgens de bovenstaande bepaling 

de genoemde a-waarden de grafiek de vorm heeft van een hoofdletter V. 

superstabiel aantrekkend dekpunt. Wanneer 3 < a < l+Vö, is er een sta- 
biele 2-cyclus. Zijn de dekpunten p en dan is de gemiddelde groei gelijk 
aan het meetkundig gemiddelde van U - 2p\ en |1 - 2^1. Een betrekkelijk 
eenvoudige berekening leidt tot het resultaat Vja^ - 2a - 4], hetgeen gra¬ 
fisch vertaald kan worden in bogen van een cirkel en van een hyperbool. 
Daarna kunnen we alleen maar experimenteel verder komen. Tot aan de 
Feigenbaum-grens van 3.669 zijn er alleen maar stabiele cycli met een 
steeds hogere periode. Maar daarna wordt het pas echt de moeite waard. 
Theoretisch verwachten we tot aan a = 4 meestal chaotisch gedrag met 



vervolgen we dat bolletje op zijn weg langs de baan, dan treedt in de af¬ 
richting wordt het volume-elementje uitgerekt terwijl in een andere rich¬ 
ting verkorting kan optreden. Gaan we niet te ver dan heeft het bolletje 
bij benadering de vorm van een drieassige ellipsoïde aangenomen. Nu 
wordt de vorm van een dergelijke ellipsoïde gekenmerkt door de (halve) 

lengte van een halve as tot die van de straal van het bolletje. D^t levert 
drie verhoudingsgetallen met de betekenis van plaatselijke groei. Er zijn 
dus in elk baanpunt drie lokale groeigetallen, even zoveel als het aantal 
dimensies van het systeem. Een middelingsproces langs een willekeurige 
baan, in dit voorbeeld de aantrekker van Lorenz, leidt daaruit drie 
Lyapunov-getallen af als meetgrootheden van pmiddelde groei. De 

king lichten we in figuur 16.5 het principe toe van een tweetal Lyapunov- 
getallen voor een tweedimensionale afbeelding. We zien hoe een cirkeltje 
vervormd wordt tot een ellipsje als een samenspel van uitrekken en in- 

Een goed voorbeeld van een tweedimensionale afbeelding is de bak¬ 
kerstransformatie. De lokale groei kunnen we bepalen door te letten op 
de deformatie van een klein cirkeltje of vierkantje. Zonder rekenen zien 
we meteen dat het elementje in de x-richting met de factor 2 wordt uitge¬ 
rekt en in de ^'-richting een inkrimping met dezelfde factor ondergaat. 
Een cirkeltje met straal r is dus een ellipsje geworden met halve assen 2r 
en r/2. De lokale groeigetallen zijn overal hetzelfde, zodat ook de gemid¬ 
delde groei met de twee Lyapimov-getallen 2 en 0.5 kan worden beschre¬ 
ven. De overeenkomstige Lyapunov-exponenten zijn log 2 en -log 2. 

Bij de kat van Amold is de berekening iets moeilijker omdat de vervor¬ 
ming van een cirkeltje een ellips in scheve stand tot gevolg heeft. Het be- 
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WINDOWC-2.-1.5)-(2.1.5) 
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DISTl=SQR((X-Xi)*(X-Xl)+(Y-Yl)*CY-Yl)) 
DIST2=SQR((X-X2)*(X-X2)+CY-Y2)*CY-Y2)) 
C1 =H2/(DIST1 *DIST1 *DIST1) 
C2=H2/(DIST2*DIST2*DIST2) 

DELU=-C 1 * (X-Xl )-C2* (X-X2) 
DELV=-C1*(Y-Y1)-C2*(Y-Y2) 



PSETCX,Y),14 

LOOP:AS=INPUT8(l) 



F=Fa)+S(D*T 

Xa)=S*Ra)^COS(F): Ya)=S*Ra)*SIN(F) 
IF T>0 THEN CIRCLE (XV(J),YVCJ)),.3,0 
IF T>0 THEN PAINT (XVa)>YV(J)),0 
CIRCLE (Xa),Y(J)),.3,COLa) 

PAINT (X(J),Ya)),COLa) 

xva)=xo): Yva)=Ya) 


PLANEET2 





Xl=RI*COSCC*T) : Y1=R2*SIN(C*T) 




XV=COS(T'H): YV=SINCr-H) 


X1V=R1*C0S(C*(T-H)): YIV=R2*SIN(C*(T-H)) 
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